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Abstract. Spatial phase-referencing in dual-field optical interferometry is reconsidered. Our analysis is based on 
the 2-sample variance of the differential phase between target and reference star. We show that averaging over 
time of the atmospheric effects depends on this 2-sample phase variance (Allan variance) rather than on the true 
variance. The proper expression for fringe smearing beyond the isoplanatic angle is derived. With simulations of 
atmospheric effects, based on a Paranal turbulence model, we show how the performances of a dual-field optical 
interferometer can be evaluated in a diagram 'separation angle' versus 'magnitude of faint object'. In this diagram, 
a domain with short exposure is found to be most useful for interferometry, with about the same magnitude limits 
in the H and K bands. With star counts from a Galaxy model, we evaluate the sky coverage for differential 
astrometry and detection of exoplanets, i.e. likelihood of faint reference stars in the vicinity of a bright target. 
With the 2mass survey, we evaluate sky coverage for phase-referencing, i.e. avaibility of a bright enough star for 
main delay tracking in the vicinity of any target direction. 
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1. Introduction 

Phase-referencing in ground-based optical interferometry 
has been proposed for faint object imaging, for narrow- 
angle astrometry and for improved precisions in visibility 
measurements ( Mozurkc wich et al. 19881 IColavita~1 992 
IQuirrenbach et al. 1994| ). Phase-referencing is the inter- 
ferometer counterpart of Adaptive Optics (AO) developed 
for single-pupils. Inhomogeneities in the atmosphere cause 
the optical path difference (OPD) between two telescopes 
to fluctuate over a range of time scales, from a few millisec- 
onds upwards, with variance rms of about 10 /im. These 
fluctuations can be measured continuously on a bright star 
and used to track a faint star in the immediate surround- 
ing. Then faint objects can be observed and small visibil- 
ity amplitude measured with sufficient Signal/Noise ratio. 
Furthermore, differential OPD fluctuations (AOPD) be- 
tween target and reference stars should average out over 
long integrations, leading to differential phase structure 
for faint objects partially resolved, or to high precision 
astrometry for unresolved sources. 



It has been noticed that AOPD residuals over long 
integrations are proportional to the measured angle 
l|Shao & Colavita 19 92). so that high precision narrow- 
angle astrometry needs to be performed with reference 
sources in the immediate surrounding of the source target, 
i.e. within the isoplanatic patch. Angular anisoplanatism 
in optical interferometry has been recently investigated 
| |Esposito et al. 2000| ). The authors define an "isopistonic" 
angle where the interferometric AOPD (rms value) re- 
mains smaller than a tenth of the observed wavelength. 
Their analysis is based on the true variance of differen- 
tial piston. A detailed analysis of the AOPD spectrum 
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has been performed by d'Arcio (1999) who considers also 
its jitter on short exposures, but does not investigate its 
potential use in phase-referencing. 

With the condition "differential phase rms smaller 
than 1 radian" , the half-cone angle of Spatial Phase- 
Referencing (SPR) is found to be the isoplanatic angle for 
AO IjFried 1982jl . divided by (2 3/5 ) due to the two entrance 
pupils of a long baseline interferometer with uncorrelated 
phase corrugations (Colavita 1992). Taking into account 
the aperture size, differential phase excursion remains 
bounded even on separation angles exceeding slightly the 
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AO isoplanatic angle, and SPR with long integrations is 
possible. On larger angular separations, random walk of 
differential phase fluctuations prevents long term aver- 
aging. The measurement of differential piston variations 
is still possible, but with short integrations, shorter than 
some coherence time for AOPD fluctuations. Sensitivity 
is reduced as compared to small angles, but the chance 
of finding a bright enough star for phase-referencing is in- 
creased. Altogether, SPR in optical interferometry has to 
be investigated in the whole telescope Field-of-View, that 
is within a few arcminutes. 

An optimization of SPR requires a detailed analysis of 
the temporal behaviour of differential phase. Our analysis 
is based on the 2-sample (or Allan) variance of time series 
of the measured differential phase. This quantity tells us 
how AOPD variations can be recovered before being av- 
eraged. Ultimately, the true AOPD variance will give the 
achievable precision on long integrations. 

Here we disentangle the problem between differential 
piston fluctuations and single-pupil wavefront distortion, 
that is between the first mode and higher modes in the 
Zernike expansion of phase corrugations. With large aper- 
tures, each pupil AO should correct for a wide field of 
view, with significant gain in the faint source direction. 

The AOPD model of our analysis is presented in the 
next Section, and quantities relevant to phase-referencing 
are recalled and clarified: the 2-sample variance and the 
attenuation factor due to fringe smearing. Results of nu- 
merical simulations are given in Sect. 3, with application 
to the 1.8m and 8m telescopes at Paranal. We further 
discuss the optimization of exposure duration and its de- 
pendence on faint star magnitude and separation angle 
(Sect. 4). Different regimes of phase-referencing are out- 
lined in Sect. 5, together with sky coverages. 



2. Model for phase-referencing 

2.1. Power spectra of atmospheric turbulence 

A detailed analysis of power spectra of quantities relevant 



to optical interferometry can be found in Conan (2000) 



where various analytical forms are also proposed for inte- 
grals of power spectra. In this Section, we recall the main 
expressions. 

The atmospheric turbulence, assumed to be stratified 
with altitude, is described by a distribution of discrete 
turbulence layers. The strength of each layer k is char- 
acterized by an optical turbulence factor Jf., sum of the 
refractive index structure parameter Cff(h) along the line 
of sight through this layer. The power spectrum of atmo- 
spheric piston over a single aperture, from a single turbu- 
lence layer, with a von Karman turbulence spectrum and 
outer scale L , is: 



W P (f) = 0.00969 (f + / 2 )- 11/6 J k F D (f) 



(1) 



where / = I/Lq, and Fd(/) is the spatial filter due to 
finite aperture size, and 



FdU) = [2Ji(7rD/)/7r£'/] 2 for a filled circular aperture 
with diameter D. 

For two sources at angular distance 0, and a turbu- 
lent layer at distance h along the line of sight, the power 
spectrum of differential piston over a single aperture is 
obtained from the covariance function: 



W A p(f) = 2Wp(f)(l - cos(27rf • R)) 



(2) 



with R = h®. 

Similarly, the power spectrum of interferometric piston or 
OPD, with a baseline B is: 



W p D (f) = 2W P (f)(l - cos(27rf • B)) 



(3) 



and the power spectrum of differential interferometric pis- 
ton (AOPD) is obtained with the product of the two filter 
functions: 

W AOP D{f) = 4W F (f)(l-cos(27rf •R))(l-cos(2vrf -B))(4) 

With the Taylor hypothesis, turbulence is displacing as a 
whole with velocity in the turbulent layer k. In each 
layer, we take the x-dircction along the wind direction 
(f • Vfc = f x Vk). The temporal power spectrum of fre- 
quency v is, for that layer: 



1 f°° v 

Wk VkJ Wk ^Vk' ^ d ^ v 



(5) 



The temporal power spectra of independent contributions 
are summed together giving the total power spectrum 
W(v) of AOPD, for the whole atmosphere crossing. 

2.2. Relevant quantities in phase-referencing 

In the following, we consider that the main OPD is tracked 
exactly on a bright star and does not contribute to inter- 
ferometric noise. The ability to follow AOPD with time, 
i.e. to measure its variation from one sample to the next, 
is a critical issue in case of large phase excursion, i.e. 
for angular separation larger than some "isopistonic" an- 
gle. Indeed, sample phase measurements have to be un- 
wrapped for long term averaging of the differential piston. 
The relevant physical quantity is not so much the magni- 
tude of differential phase (due to differential piston), but 
each sample departure from a local average. The most use- 
ful quantity for describing such a process is the 2-sample 
variance, or Allan variance used for the characterisation 
of frequency standards (Rutm an 1978|l . The local average 
is simply made with two successive samples. Let To be the 
exposure duration for each sample AOPD^, and also the 
sampling interval. The 2-sample variance of the time se- 
ries is <j\(t ) = \ ((AOPD J+ i - AOPD,) 2 ), and it can 
be expressed in terms of the true variance Iaopd( t ) by 
ilRutman 1978jl : 



cr i( T o) = ^AOPD^o) - IaOPd( 2t o)} 



with: 



'AOPD 



(r) 



W{y)[ 



sin(7wr) , , 
-]dv 

TTVT 



(6) 



(7) 
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A third relevant quantity is the phase jitter within 
each sample, responsible for fringe smearing and then sig- 
nal loss. A detailed analysis of fringe smearing in VLBI, 
due to instabilities in frequency oscillators, can be found 



Thompson et al. (2001) A similar approach applied 



to differential piston in dual-field optical interferometry 
yields the following expression for the attenuation factor 
due to finite exposure duration. 

Let A<fi(t) be the differential interferometric phase (due 
to AOPD), and to the exposure duration. The (random) 
attenuation factor gi and measured differential phase A<f>i 
are given by: 



die 



1 

TO 



Li+ra 



(8) 



With a quadrature-phase scheme for fringe phase and am- 
plitude measurements, the average attenuation factor use- 
ful for Signal/Noise ratio estimate, G(to), is taken as the 
square root of the quadratic average of gi : 



G(r ) = 



1 



'o Jo 



(exp j[A<t>(t) - A(j){t')])dtdt' (9) 



1/2 



If A(j> has gaussian statistics and is a stationary process, 
the quantity between brackets is expressed in terms of 
the structure function of differential interferometric phase: 
°A<*( r ) = 2 (x) 2 [ i? (°) - i? ( r )]> where A is the optical 



wavelength, and R(t) the correlation function or Fourier 
Transform of W(v): 



R(t) = 2 W{v) 
J o 



cos(27r vr)dv 



Equ. 9 is now written: 

G(t ) = — (1 )exp( ^ 

Lto Jo t o 2 



%dt} 1/2 



(10) 

(11) 
for 



A similar expression is given by Buscher 
the smearing factor on main OPD, the fringe corre- 

(T 2 (t) 

lation function exp( — ^| — ) being expressed in terms 
of an interferometer coherence time. In previous works 
on off-axis phase-referencing, the fringe smearing factor 
has been often approximated by exp(— 27t 2 cr^ es (r ) /A 2 ), 
where af es (ro) = 2 W(u)[l — sinc(7T^r)] 2 ciz/, valid only 
in the weak smearing limit, that is useless for phase- 
referencing beyond the "isopistonic" angle. 

2.3. Turbulence model of the atmosphere 

Distribution profile of turbulence parameters above 
Paranal are taken from the results of the PARSCA cam- 
paign of March 1992 l|Fuchs & Vernin 1993jl . with bal- 
loon soundings and SCIDAR measurements. The aver- 
age turbulence profile [C 2 (/i)] and wind amplitude pro- 
file [V(/i)] have been used by several authors for sim- 
ulations of atmospheric effects < |Delplancke et al~ 000 
Esp osito et al. 2000| ), with a nominal seeing of 0.65". 
Continuous profiles from the PARSCA report have been 
re-sampled with 1 km step (which is about the vertical 



resolution of SCIDAR measurements). Turbulence content 
summed over each layer is the optical turbulence factor Jk 
used in the following. Wind velocity has also been aver- 
aged over each layer and is shown on Fig. 1. In this model, 
the weighted average wind velocity equals 8.5m/s. 




5 10 
Altitude above Paranal [km] 



Fig. 1. Distribution of turbulence parameters with alti- 
tude above Paranal, from the average model of PARSCA 
1992 campaign. The turbulence factor is in m 1 / 3 . The 15 
layers adopted in our model are marked with triangles. 



Wind direction is another useful parameter but is not 
available in the report on the PARSCA campaign. Wind 
direction changes somehow with altitude, and scatter in 
wind direction smoothes singularities of the "frozen in" 
turbulence model. We have taken a dominant wind direc- 
tion mainly perpendicular to the interferometer baseline, 
that is a situation where turbulence effects in long baseline 
interferometry are the most significant. 

Another parameter of the model is the outer scale of 
turbulence in each one of the layers (Lo,fc). This quan- 
tity is hardly measured locally. The Generalized Seeing 
Monitor (GSM) can measure a weighted average of Lo,fe 
along the line of sight. Strong turbulence layers with small 
Lq values have a dominant effect in this weighted average 
HBorgnino~1 990 ) . The spatial coherence outer scale mea- 
sured in this way at Paranal has a log-normal distribution 
with a median value of 22m l|Martin et al. 2000|l which 
may reflect mainly turbulence properties close to the tele- 
scope. The outer scale high in the atmosphere has a signif- 
icant impact on long term averaging, and a limited impact 
on short term averaging of differential phase, as shown in 
the next Section. So, the lack of a detailed knowledge of 
the outer scale distribution with altitude does not impinge 
on our study. 
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; L =250m ; 



10 r main OPD 100m 




Exposure duration [sec] 

Fig. 2. Main OPD variance, and AOPD variance for sep- 
aration angles 15" , 30" and 60" , 100m interferometer 
baseline and 8m aperture diameter. True variance and 
2-sample variance are respectively shown with thin and 
thick lines. The other parameter is the outer scale of tur- 
bulence: 50m (dashed lines), 100m (full lines) and 250m 
(dash-point lines). For 30" angular separation, only the 
full line is shown for clarity of the Figure. 

3. Averaging of atmospheric effects: impacts of 
separation angle and outer scale 

Simulations of atmospheric turbulence have been per- 
formed with an interferometer baseline of 100m and tele- 
scope diameters of 8m and 1.8m. In our model, the sepa- 
ration vector between the two stars is parallel to the base- 
line direction, and the zenithal distance is nearly zero. 
The power spectral densities (main OPD and AOPD) for 
each one of the 15 turbulent layers of Fig. 1 are estimated 
with Equ. 3-5, then the true variance and 2-sample vari- 
ance, summed on the whole line of sight, are estimated 
with Equ. 6-7 for different exposure durations r. Results 
are shown in Fig. both for main OPD and for AOPD. 
Impacts of separation angles 9 ranging from 15 to 60" 
and of outer scale Lq between 50 and 250 m are shown. 
For short exposures (r < a few 0.1 s), differential phase 
variations, from one sample to the next, is small enough at 
optical wavelength to allow reconstruction of a continuous 
and unambiguous differential phase solution from succes- 
sive samples. The relevant quantity for differential phase 
reconstruction is then the 2-sample variance, and not the 
true variance. This difference is essential since the true 
variance is a monotonic decreasing function with expo- 



sure duration whereas the two-sample variance increases 
first, then peaks at r max ~ 1-2 seconds, and ultimately 
approaches the true variance for long exposures. 

In Fig. 2, the variance is found to be proportional 
to the square of the angular separation, as already 
shown with analytic expressions (Sha o fc Colavita 1992 
IColav ita 1994). We point also that: 

— the outer scale of turbulence has a significant impact 
both on main OPD, and on AOPD for long exposures, 
and a small impact on the 2-sample variance for short 
exposures, 

— the variance of AOPD remains smaller than 1/im 
rms, or much smaller than the coherence length of 
starlight in the near-IR (with less than 10% relative 
bandwidth). This makes it possible differential delay 
tracking without feedback, that is open loop operation 
whatever the separation angle and exposure duration. 



Table 1. Greenwood time delay ^o, isoplanatic angle 6>o, 
and "isopistonic" angle 9\ for 1.8m (AT) and 8m (UT) 
telescopes in the H and K band, at zenith, at Paranal. 



IR band 


H 


K 


t [ms] 


20 


29 


00 


8.6" 


12.1" 


01, AT 


10.7" 


14.3" 


01,UT 


17.4" 


23.0" 



The Greenwood (or AO) time delay is to = 0.314 ro/V 
(Frie d 1990jl . where r is the Fried parameter and V the 
average wind velocity. For single aperture, the isoplanatic 
angle 6 = [2.91(f) 2 £ fc J k h]/ Z ]~^ ^Frieri 198211 is de- 
fined as the angular extent over which differential wave- 
front phase corrugations is smaller than 1 radian (rms). 
By analogy, in interferometry, an angle 9\ can be defined 
as the angular range where 27rAOPD/A (rms value of the 
true variance) is smaller than 1 radian. The estimation 
of this parameter is in Table \I\ for two different values of 
telescope diameter and for an outer scale of 100m. The 9q 
and 9i : ut estimates made by Esposito et al. (2000) for the 
same observing site are somewhat smaller due to a more 
pessimistic turbulence model and to a more conservative 
definition of the "isopistonic" angle (rms phase of 27r/10 
instead of 1 rad.). 

The average smearing factor G(tq) given by Equ. 11 
is shown in Fig. for AOPD measurements with 4 sep- 
aration angles, from 15 to 60". This smearing factor in 
phase-referencing remains close to 1.0 for sample exposure 
shorter than 0.1s, for the whole range of angular separa- 
tion. For large enough angles and long exposure duration 
t, we verify that the smearing factor is in t -1 / 2 as shown 
by Buscher lfT58"8"jl . 
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before detection, the source noise contribution is further 
reduced and we get: 



Fig. 3. Smearing factor due to finite exposure duration 
for differential OPD measurements with separation angles 
15" , 30" , 45" and 60" , in the H and K bands (respectively 
dashed and full lines). The smearing factor is also shown 
for main OPD measurements. 

4. Sensitivity estimates in phase-referencing 

Differential piston measurements is affected by the limited 
Signal/Noise ratio that can be achieved on faint sources. 
With B and S respectively the noise and signal in visi- 
bility measurement of the faint source, the measured 2- 
sample differential OPD variance is the sum of two con- 
tributions, a measurement noise and an atmospheric dif- 
ferential OPD: 



Y 2 (t ) = a 2 m (r ) + a 2 A (r ) 



The measurement noise variance, given by: 

A, 2r B(r ), 
2ir } 



~* In- \ I " \ 2 / "VV \2 



(12) 



(13) 



is a decreasing function of the coherent exposure duration 
To, whereas the 2-sample variance a\ of AOPD first in- 
creases from zero to its peak value, and then decreases. 
The quadratic sum of the two often shows a local mini- 
mum that gives an optimum exposure duration. 

4.1. Signal and Noise 

Sensitivity estimate is performed in the H and K bands. 
Measurement noise contributions are assumed to come 
from read-out noise of the detector, and from star and 
background photon noise. Let oi be the read-out noise, 
in electrons per measurement, 60 and sm the number of 
photons that would be detected per second and per aper- 
ture from sky background and from a star with magni- 
tude M, respectively, without mode filtering. Fringe phase 
measurement is supposed to be performed with an ABCD 
scheme and 4 simultaneous detector readouts, also known 
as 'spatial discrete modulation scheme with 4 phases' 
< |Cassaing et al. 2000| ). In the photon-rich limit, the Signal 
and Noise of such a detection scheme have been derived 
( Dai gne fc Lestrade 1999| ) . With additional mode filtering 



B 2 (t ) = 2<j? + (bo + Ts M )T 



(14) 
(15) 



7 being the interferometer visibility, T the average Strehl 
ratio or aperture efficiency in the point source direction, 
and G(tq) the average fringe smearing factor for expo- 
sure duration To (Equ. 11). Strehl ratio fluctuations, 
uncorrelated on the two pupils, are responsible for an ad- 
ditional loss of coherence due to unbalanced amplitudes, 
which is not considered here. 

Parameters of SNR estimates for the UTs are given 
in Table |21 The overall photometric efficiency is the arm 
transmission factor (0.30 and 0.31 in the H and K band 
respectively) times additional factors: coude train with 
two-field selection and AO, supposed to be 0.95 in both 
bands, overall efficiency of FSU (optics, fiber injection and 
mode filtering, detector quantum efficiency), supposed to 
be 0.3 in each band (Cassaing, private communication). 
The Strehl ratio is taken from the MACAO document, 
and the interferometer visibility is an estimate. The read- 
out noise is taken as 7 electrons per measurements, with 
4 frequency channels per band (spectral resolution ~35). 



Table 2. Parameters for Signal and Noise estimates with 
UTs. _Ei5 is the reference photon flux for a source of magni- 
tude 15 ffrom lCohen et al. 19921 at Mauna Kea). S15 is the 
product of £15 by the throughput (collecting area times 
bandwidth, times overall photometric efficiency). The col- 
lecting area of each 8m telescope is 49.1 m 2 . 



IR band 


H 


K 


wavelength [/xm] 

E15 photons/ (m 2 .s.^m) 


1.65 


2.2 


9560 


4560 


bandwidth \pm] 


0.19 


0.25 


overall photometric efficiency 


0.085 


0.088 


sis [photo-e~/s] 


7581 


4925 


sky background [mag./C')^] 
beam size [(") 2 ] 


14.4 
0.01 


13.0 
0.01 


60 [photo-e~/s] 


132 


310 


r Strehl ratio 


0.30 


0.45 


7 Interf. visibility 


0.9 


0.9 


S15 coherent [photo-e~/s] 


2047 


1995 



A similar SNR estimate is obtained with a single fre- 
quency channel per band, and a read-out noise of 10 elec- 
trons instead of 7. Spectral resolution (e.g. for group delay 
measurements with a single band) is not so much a penalty 
in terms of sensitivity threshold. Despite the smaller Strehl 
ratio at H band, the coherent signal from an unresolved 
star will be larger than at K band. With fewer photons 
from sky background, phase-referencing at H band is most 
useful, as shown in the next section. 
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Fig. 4. The 2-sample phase rms of main OPD and stars 
with magnitude 11-15 and 13-17 versus exposure duration 
with VLTI auxiliary telescopes (ATs) and main telescopes 
(UTs), in the H and K bands (respectively dashed and 
full lines). The outer scale of turbulence is Lo=100m. The 
sensivity gain with UTs is about 2 magnitudes compared 
to ATs. 

Sensitivity estimates have been performed also with 
ATs. Starting from the UT's parameters (Table 2), col- 
lecting area is divided by 20, beam size is multiplied by 
10, and Strehl ratios are taken as 0.54 and 0.72 in the H 
and K bands respectively (Tip-Tilt corrections only). The 
other parameters (spectral resolution, read-out noise . . . ) 
are kept unchanged. 




10 1 10 

Exposure duration [sec] 

Fig. 5. The 2-sample phase rms of main OPD (thick 
lines) and AOPD (thin lines) in the H and K bands (re- 
spectively dashed and full lines), with 8m telescopes and 
Lq = 100m, for three values of separation angle 9. On each 
graph, the magnitude range for AOPD is from 15 to 19. 
A single magnitude is shown for main OPD. 



slightly decreases with separation angle. The range we find 
for an optimum exposure is from 0.03 to 0.5 seconds. 



4.2. Magnitude limits 

The 2-sample phase rms reached with piston and measure- 
ment noise is: 

y(r Q ) = ^Y(to) (16) 

Main OPD phase is first plotted in Fig. for star 
magnitudes ranging from 11 to 17. A minimum rms value 
is obtained for an exposure duration r m . The exposure T op t 
which optimizes the SNR per time unit will be slightly 
smaller than r TO . The minimum phase rms is about the 
same in the two bands, and stars as faint as H and/or 
K = 15 might be used for main delay tracking with 8m 
telescopes. Close loop operation will degrade such a figure 
of merit, so that the magnitude limit will most probably 
lie in the range 13-14. 

The 2-samplc differential phase rms measured with the 
UTs is plotted in Fig. [jj for separation angles 15" , 30" and 
60", and for faint star magnitude in the range 15 to 19. 
Again a minimum is reached for stars not too faint. The 
exposure time at minimum increases with magnitude, and 



5. Phase-referencing with short exposure 

The minimum phase variation reached for r = r TO is about 
the same in the H and K bands, so that we should consider 
the benefit of two-band observing. Let Co be the minimum 
phase variation (rms) at the effective wavelengths \h and 
\k- At minimum, the contributions of differential piston 
and measurement noise to the measured 2-sample variance 
are about the same (Equ. 12). Measurement noises in the 
two bands being uncorrelated, the uncertainty on AOPD 
variation, from one sample to the next, is: 

v sl ~ + A|-) 1/2 (17) 

For <7o = l radian, the uncertainty on AOPD variation is 
about 0.21 /im that is less than A/10 in the K band. Such 
a figure is obtained with an equivalent SNR of about 1.6. 
Unwrapping differential phase measurements, from one 
sample to the next, should be possible with filtering al- 
gorithms reducing the contribution of measurement noise. 
Indeed, in the considered frequency range of about 10 Hz, 
the AOPD power spectrum will have a steep slope due to 
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Fig. 6. Different regimes of Phase-Referencing based on 
the 2-sample differential phase. In the (a) domain, there 
exists a well defined minimum satisfying AOPD phase 
rms smaller than 1 rad for an optimum exposure dura- 
tion T op t < Is. Domain (b) corresponds to short expo- 
sure without a well defined optimum exposure duration. 
Domain (c) corresponds to long integration, limited in the 
upper right corner by a straight line with a somewhat ar- 
bitrary threshold on phase rms, taken here to be smaller 
than 1 rad for 20 seconds integration. 

spatial filtering with the telescope aperture, whereas mea- 
surement noise will have a flat spectrum. Its contribution 
could be reduced with some time averaging. 

A domain of SPR with short exposure is defined in 
a diagram (separation angle, magnitude of the faint ob- 
ject), with the condition: "there exists a minimum in the 
2-sample AOPD phase, and the minimum value is smaller 
than 1 radian (rms)" . An optimum exposure duration for 
faint source observing is then clearly defined. This domain 
is shown with label (a) in the lower part of Fig. H3 For 
fainter sources and separation angle not too large (< 30"), 
there may be no minimum, but a measured 2-sample phase 
rms smaller than 1 radian for exposures smaller than 
about 1 second, as, for example, K > 18 with 9 — 15" 
(Fig. 5). Phase-referencing is still possible with short ex- 
posure, but without an optimum duration. Uncertainty is 
then dominated by measurement noise, and this domain is 
labelled (b) in Fig.[fj] The third domain, labelled (c), is the 
usual domain for long exposures. The exposure duration 
should be long enough in order to smooth out fluctuations, 
that is well beyond the peak of the 2-sample variance. As 
shown in Fig. [5] phase-referencing of a 19th magnitude 
source is marginally possible up to a separation angle of 
1 arcminute, the phase rms being about 1 radian for 20 
seconds integration. 

The newly identified domain (a) for SPR with short 
exposure is useful only if there are stars satisfying the 
magnitude and separation constraints shown in Fig. We 
first consider phase-referencing for the astrometric search 
of exoplanets orbiting bright stars, i.e. bright enough for 
main delay tracking. We would like to measure the angu- 
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Fig. 7. Accessibility to faint stars in terms of magnitude 
and separation angle in the H and K bands (respectively 
dashed and full lines). The 3 domains of Fig. 6 for SPR 
are recalled here with thick lines for ATs and UTs obser- 
vations. The upper right limit of the domain (c) is omitted 
for clarity. The thin curves indicate that there are 3 stars 
in average within the half cone of the corresponding sep- 
aration and brighter than the corresponding magnitude. 
These curves are labelled with the galactic latitude of the 
target direction. 

lar separation between a main target and reference stars in 
its surrounding. As randomly choosen stars can turn out 
to be binaries, we take a conservative value of 3 reference 
stars (in average) within the useful solid angle. Reference 
stars for differential astrometry will most often be faint 
objects, so that we must rely on star count models for 
their distribution with galactic coordinates. Our results 
are based on the Besangon model on stellar populations 
m our Galaxy jReyle fc Robin 2001) . The half cone angle 
with 3 stars (in average) brighter than a given magnitude 
m is taken as [7V(m)] _1 / 2 , where N(m) is the cumulative 
density of stars brighter than m. Results shown in Fig. [7] 
are for a region of the sky suited for phase measurements, 
with a nearly zenith transit at Paranal. Its galactic longi- 
tude is —15 degrees, and the latitude range is from —15 
to —60 degrees. Whatever the target direction, 3 stars are 
likely to be found with UTs and SPR with short expo- 
sure (domain (a)). At a galactic latitude of —60 degrees, 
the separation angle reaches 1 arcmin and faint star mag- 
nitude 17.5. Long exposure SPR with fainter and closer 
stars is also possible, say 19.5 magnitude and about 42", 
but in a regime where the useful photon flux from the star 
is one tenth the photon flux from the sky background in 
the K band. Closer to the galactic plane, SPR with short 
exposure is much easier due to the greater star density. 

A different estimate has to be performed for phase- 
referenced imaging of a faint object. What is needed is 
a bright star in its surrounding for main delay tracking. 
The bright (or guide) star lies in the magnitude range of 
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Fig. 8. Probability of finding a guide star for main delay 
tracking in the vicinity of any target direction with galac- 
tic latitudes 15, 30, 45 and 60 degrees. The guide star has 
to be brighter than K=11.0 (top), 13.0 (middle) or 14.0 
(bottom). 



near-IR surveys, and our estimate of sky-coverage is based 
on the Second Release of the 2MASS survey 1 . The prob- 
ability of finding at least one star brighter than Ku m has 
been estimated in terms of the size of the solid angle for 
phase-referencing, at different galactic latitudes (Fig. |SJ| . 
Results on positive and negative galactic latitudes have 
been averaged, and the magnitude limit for the guide star 
is Ku m = 11 (with two ATs), and 13 or 14 (with two UTs). 
Taking a galactic latitude of 30 degrees as an average sky 
density, the probability of finding a guide star with 2 ATs 
(Kum < 11, within 20") is about 2%. Significant sky cov- 
erage, with a probability of 50%, will be reached only with 
2 UTs and improved corrections for the atmospheric tur- 
bulence {Kum < 14 and AO corrections in two different 
directions, up to 1' apart). 



1 This publication makes use of data products from the the 
Two Micron All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and the 
National Science Foundation. 



6. Discussion 

The useful angular range for phase-referencing and dif- 
ferential astrometry is not strictly limited, but dependent 
on the magnitude of the faint object (Fig. 0). Then, the 
'isopistonic' angle that has been defined in the literature is 
too restrictive. In our analysis, there exists a specific angle 
given by the 'triple-point' in Fig. |BJ where the 2-sample 
phase rms equals 1 radian at a local minimum. The sep- 
aration angle and the magnitude of the faint object at 
this triple-point depend both on the interferometer sen- 
tivity and on the site turbulence. It characterizes a whole 
system (atmosphere + telescope + instrument). An inter- 
esting finding of our calculation is that the triple-point is 
located at about (20" , 16) for ATs and (30" , 18.3) for UTs 
with similar instruments. This makes a very significant dif- 
ference in sky coverage, differential astrometry with ATs 
being limited to rather dense stellar fields (Fig. [TJ) . 

Short exposure measurements of differential phase 
fluctuations integrate various contributions: differen- 
tial piston, sky background variations, interferome- 
ter/instrument instabilities ...Short term measurements 
can then be edited, filtered and weighted before being av- 
eraged to properly extract source structure and relative 
position information. The similar sensitivities and star 
counts obtained in the H and K bands (Fig. 0) must be 
stressed, as dual-band observations allow more robust pro- 
cedures than single band near the magnitude limit. 

In our analysis, Adaptive Optics has not been ad- 
dressed, whereas wavefront corrections are needed in 
the two sources direction, at least with large telescopes. 
Tomographic sensing of the atmospheric turbulence with 
the observation of several guide stars has been pro- 
posed for wavefront corrections on a whole telescope FOV 
dTallon fc Foy 1990D . The technique is known as Multi 
Conjugate Adaptive Optics (MCAO). Although mainly 
concerned with improved Strehl Ratio, MCAO is supposed 
to correct for atmospheric turbulence in a whole volume, 
and hence to correct for differential pathlength fluctua- 
tions on any entrance pupil. Differential interferometric 
piston (AOPD) should be significantly reduced, allowing 
for long term averaging in faint object imaging on a whole 
telescope FOV. However, it is difficult to know how much 
AOPD would be reduced since corrugations with large 
scale lengths, about the pupil size, are less efficiently cor- 
rected than shorter scale lengths. 

Sky coverage with natural guide stars was thought to 
be a main problem for AO and not too large telescopes. An 
efficient technique for MCAO has been proposed recently 
agazzoni et al. 2002| ). It is layer-oriented, with multiple 
Field of View: a large annular FOV is used for sensing low 
altitude turbulence, whereas a narrower FOV (about 2') 
in the telescope axis is used for sensing high altitude tur- 
bulence. Sky coverage exceeding 50% should be reached 
with 8m telescopes and a wavefront sensor in the R band. 
Such a sky coverage is quite consistent with our estimate 
on spatial phase-referencing in dual-field optical interfer- 
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ometry in the near-IR. It makes faint source imaging with 
VLTI/UT a very promissing next development step. 

With the AO system being presently implemented at 
the Coude focus of UTs (MACAO), faint source imaging 
will be restricted to the isoplanatic patch of guide stars. It 
may reach the telescope FOV of 2 arcmin when turbulence 
in the upper layers of the atmosphere is particularly weak. 
For bright enough targets, our analysis shows that AOPD 
will be measured and tracked well beyond the isoplanatic 
angle. 
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